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(57) Abstract 

The invention relates 
to a method for detecting 
modes of a dynamic system 
comprising a plurality of 
modes having respectively 
a set o-(t) of characteristic 
system parameters. According 
to this method, a chronological 
series of at least one system 
variable x(t) is subject to a 
modelling, such as a switch 
segmentation, designed 
to detect a predetermined 
prediction model, such as 
a neuronal network relative 
to a corresponding system 
mode, in each time interval 
with a predetermined minimal 
length for each system variable 
x(t). After modelling the 
chronological series, a drift 

segmentation is performed wherein, a sequence of mixed prediction models is detected in each time interval where the system shifts from 
a first system mode to a second system mode, said sequence being a result of a linear twin superposition of the prediction models of both 
system modes. 
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(57) Zusammenfassung 

Bei einem Verfahren zur Erfassung der Moden eines dynamischen Systems mit einer Vielzahl von Moden, die jeweils einen Satz <*(t) 
charakteristischer Systemparameter besitzen, wird eine Zeitreihe mindestens einer Systemvariablen x(t) einer Modellierung wie z. B. einer 
Schaltsegmentierung unterzogen, die dazu eingerichtet ist, in jedem Zeitabschnitt einer vorbestimrnten Mindestlange fiir jede Systemvariable 
x(t) ein vorbestimmtes Vorhersagemodell wie z. B. ein neuronales Netzwerk fiir eine entsprechende Systemmode zu erfassen, wobei nach der 
Modellierung der Zeitreihe eine Driftsegmentierung erfolgt, bei der in jedem Zeitabschnitt, in dem das System von einer ersten Systemmode 
zu einer zweiten Systemmode iibergeht, eine Folge von gemischten Vorhersagemodellen erfaBt wird, die durch eine Iineare, paarweise 
Oberlagerung der Vorhersagemodelle der zwei Systemmoden gegeben ist. 
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Verfahren zur Erfassung zeitabhangiger Moden 
dynamischer Systeme 

Beschreibung 

Die Erfindung betrifft ein Verfahren zur Erfassung dynamischer 
Systeme, die durch im Zeitverlauf nichtstationare 
Systemparameter charakterisierbar sind, insbesondere ein 
Verfahren zur Segmentierung von Zeitreihen von Mefigrofien 
(Variablen) dynamischer Systeme und zur Identif izierung der 
die Segmente charakterisierenden Systemparameter (Moden) . 

Unter einem dynamischen System wird hier insbesondere jede 
Erscheinung verstanden, deren Zeitverlauf durch eine diskrete 
Abb il dung vom Typ 

x(t+l) = f a(t) (a(t) ) (0.1) 

darstellbar ist. Es werden jedoch auch Systeme mit mehreren 
(z. B. zwei) simultan auf genommenen Zeitreihen x, y gemaii 

y (t+x) = f a(tJ (x(t) ) (0.2) 

betrachtet. Dabei bezeichnen a(t) einen Satz 

charakteristischer Systemparameter, x einen Zustand, der im 
allgemeinen einen Vektor in einem mehrdimensionalen 
Zustandsraum bildet, und y einen zeitlich verschobenen 
Zustand. Der Zustandsraum wird von Variablen aufgespannt, die 
z. B. physikalische, chemische, biologische, medizinische, 
geologische, geometrische, numerische und/oder 
prozelitechnische Grofien sein konnen . 

Die Anzahl der Systemvariablen, die zusammen mit der Dynamik f 
das System beschreiben, entspricht der Dimension des 
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Zustandsraumes. Hier werden Systeme betrachtet, deren 
Parameter a ebenfalls zeitlich veranderlich sein konnen. Ein 
gegebenes System mit zeitlich unveranderlichen Parametern a 
wird im folgenden auch als Mode bezeichnet. 

Beobachtbare oder mefibare Sys temvariablen (Mefigroflen) bilden 
erfaftbare Zeitreihen oder Datenstrome, die fur die jeweilige 
Abfolge von Systemmoden charakteristisch sind. Sind die 
Systemparameter innerhalb der Zeitreihen liber bestimmte 
Zeitabschnitte (Segmente) unveranderlich, so kann die 
Zeitreihe entsprechend den jeweils bestehenden Systemmoden 
unterteilt werden (Segmentierung) und jedes Segment einer 
Systemmode zugeordnet werden ( Identif izierung) . 

Viele Vorgange in der Natur wie auch bei technischen 
Anwendungen wtirden sich dann vorhersagen und/oder 
kontrollieren lassen, wenn die ihnen zugrundeliegenden 
dynamischen Prozesse mathematisch modelliert werden konnen. 
Die Analyse und Charakterisierung praktisch gegebener 
dynamischer Systeme wird oft dadurch erschwert, da£ sich die 
Systemmoden wahrend der Beobachtung andern . Beispiele hierfiir 
sind allmahiiche Veranderungen, die sich in Drifts oder Trends 
der Systemparameter widerspiegeln, oder auch spontane oder 
abrupte Veranderungen in der Dynamik komplexer Systeme, etwa 
wenn sich Konf igurationen spontan oder von aufien getrieben 
plotzlich andern. 

Ein Beispiel fur ein betrachtetes System ist die Erzeugung von 
Sprachsignalen im Mund-Rachen-Raum, bei der das System standig 
seine Konf iguration und somit seine Mode andert. Es besteht 
ein starkes Interesse daran, die Moden zu erfassen und zu 
identif izieren, die einem beobachteten Zeitverlauf einer 
Variablen (im Beispiel: Luf tdruckschwankungen) zugrundeliegen, 
auch urn verbesserte Vorhersagen uber das betrachtete System 
treffen oder dieses besser kontrollieren zu konnen. 
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Prinzipiell konnen dynamische Systeme anhand vora gemessenen 
Signalen analysiert werden und es sind Reihe von Methoden 
bekannt, aus Zeitreihen Modelle zu gewinnen, die sich fiir die 
Vorhersage und eine Kontrolle des Systemverhaltens eignen. So 
ist bekannt, dafi sich der Zustand eines dynamischen Systems 
durch Erfassung der Zeitabhangigkeit beobachtbarer Mefigrofien 
modellieren lafit. Diese Modellierung erfolgt gemafi einem 
ersten Ansatz durch eine Rekonstruktion des Zustandsraumes 
mittels sogenannter Zeit-Verzogerungs-Koordinaten, wie es z. 
B. von N. H. Packard et al . in "Physical Review Letters" (Bd. 
45, 1980, S. 712 ff.) beschrieben ist. Auf der Basis einer 
solchen Rekonstruktion lafit sich dann nur ein einziges 

(globales) Modell f fur die Dynamik finden. Die globale 
Rekonstruktion des Systems ist ferner nachteilig, da bei 
Anwendungen fur vieldimensionale Systeme eine Vielzahl von 
Eingangsgrbfien als Randbedingungen vorab bekannt sein mussen 
und/oder aufgrund der hohen Dimensionality das System sich 
praktisch nicht mehr schatzen (erfassen, abbilden) lafit 
und/oder ein ubermafiig hoher, nicht praktikabler Rechenauf wand 

entsteht . 

Aufierdem ist dieses Verfahren fiir den Fall zeitlich 
veranderlicher Parameter allgemein nicht anwendbar . Die 
Analyse und Modellierung von dynamischen Signalen wird jedoch 
haufig dadurch erschwert, dafi die zugrundeliegenden Systeme 
sich in der Zeit in wesentlichen Parametern verandern. 
Beispiele sind z. B. Signale aus der Medizin, bei denen ein 
Organ wie das Herz Oder das Gehirn viele dynamische Moden hat, 
die sich abwechseln, oder etwa den Sprachsignalen, bei denen 
das erzeugende System, der Mund-Rachen-Raum, in der Zeit 
of f ensichtlich verschiedene Konf igurationen durchlauft. 

Gemafi einem weiteren Ansatz ist daher aus der Publikation von 
K. Pawelzik, J. Kohlmorgen und K.-R. Muller in "Neural 
Computation" (Bd. 8, 1996, S. 340 ff.) bekannt, Datenstrome 
entsprechend zunachst unbekannten, zeitlich wechselnden 



WO 99/14686 



PCT/EP98/05793 



Systemmoden durch Simuiierung mit mehreren, untereinander in 
Konkurrenz stehenden Modellen zu segmentieren. Die Modelle 
werden bevorzugt durch neuronale Netzwerke gebildet, die 
jeweils fur eine Dynamik charakteristisch sind und nach 
vorbestimmten Trainingsregeln urn die Beschreibung der 
einzelnen Punkte des Datenstroms konkurrieren . 

Mit dieser Methode gelingt es, eine Zeitreihe in Abschnitte 
quasistationarer Dynamik zu zerlegen und simultan Modelle fur 
diese System-Moden aus den Zeitreihen zu identif izieren. 

Die Segmentierung gemaft K. Pawelzik et al . , zu der unten 
Einzelheiten angegeben werden, erlaubt die Zuordnung von 
Segmenten zu bestimmten Systemdynamiken oder -moden und fuhrt 
zu einer Erfassung des Datenstroms als Vorgang mit einem 
diskreten Schalten (sog. "switching") zwischen den Moden. 
Diese Beschreibung der Parameterdynamik komplexer Systeme 
stellt zwar gegenuber der o. a. globalen Modellierung 
hinsichtlich der Genauigkeit und der Segmentierung 
verschiedener Systemzustande einen Fortschritt dar . Allerdings 
kann der Obergang zwischen verschiedenen Systemzustanden nicht 
ausreichend beschrieben werden. Es hat sich insbesondere bei 
der Analyse realer Systeme z. B. bei medizinischen Anwendungen 
gezeigt, daft die Segmentierung auf bestimmte Falle mit 
moglichst klaren Modenunterschieden und geringem Rauschen 
beschrankt, allgemein jedoch bei zeitlichen Veranderungen der 
erzeugenden Systeme unzuverlassig ist. 

Solche zeitlichen Veranderungen der erzeugenden Systeme machen 
die beobachtbaren Signale nichtstationar und ftihren dazu, dafi 
sich die Systeme im allgemeinen nicht mehr durch einheitliche 
Modelle beschreiben lassen. Erfolgen diese Veranderungen der 
Systeme plotzlich, so spricht man von Sprungprozessen . 

Die Aufgabe der Erfindung ist es, verbesserte Verfahren zur 
Erfassung der Moden dynamischer Systeme mit nichtstationaren 
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Systemparametern anzugeben, mit denen die Beschrankungen 
herkommlicher Verfahren uberwunden werden konnen und die es 
insbesondere ermoglichen, automatisch mit praktikablem 
Bearbeitungsauf wand und hoher Zuverlassigkeit Segmentierungen 
und Identif izierungen von Zeitreihen mit einer erhohten Zahl 
von Einzelheiten vorzunehmen . 

Diese Aufgabe wird durch das Verfahren mit den Merkmalen gemaft 
den Patentanspruch 1 gelost. Vorteilhafte Ausfiihrungsf ormen 
der Erfindung ergeben sich aus den abhangigen Ansprlichen. 

Die Erfindung basiert auf der Idee, Ubergange zwischen 
verschiedenen Moden eines dynamischen Systems als 
Zwischenmoden des Systems aufzufassen, die paarweise lineare 
Interpolationen der Ausgangs- bzw. Endmoden des Ubergangs 
darstellen. Die betrachteten dynamischen Systeme gehen eher 
allmahlich von einer Mode in die andere liber, statt abrupt 
zwischen Moden zu schalten. Die Erfindung zielt darauf ab, 
solche Ubergange zwischen dynamischen Moden in Signalen und 
die Moden zu identif izieren . 

Es wird daher bei einem Verfahren zur Erfassung der Moden 
dynamischer Systeme z. B. nach einer Schaltsegmentierung einer 
Zeitreihe mindestens einer der Systemvariablen x(t) des 
Systems eine Drif tsegmentierung vorgenommen, bei der in jedem 
Zeitabschnitt, in dem das System von einer ersten Systemmode s^ 
zu einer zweiten Systemmode s 3 iibergeht, eine Folge von 
gemischten Vorhersagemodellen g i erfafit wird, die durch eine 
lineare, paarweise Uberlagerung der Vorhersagemodelle f i#j der 
zwei Systemmoden s l#j gegeben ist. 

Gegenstand der Erfindung ist auch eine Vorrichtung zur 
Erfassung eines dynamischen Systems mit einer Vielzahl von 
Moden s i mit jeweils charakteristischen Systemparametern a(t). 
Die Vorrichtung enthalt eine Einrichtung zur Aufnahme einer 
Zeitreihe mindestens einer der Systemvariablen x(t) des 
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Systems, eine Schaltsegmentierungseinrichtung, die dazu 
eingerichtet ist, in jedem Zeitabschnitt einer vorbestimmten 
Mindestlange fur die Systemvariablen x(t) ein vorbestimmtes 
Vorhersagemodell f i fur eine entsprechende Systemmode s i zu 
erfassen, und eine Drif tsegmentierungseinrichtung, mit der in 
jedem Zeitabschnitt, in dem das System von einer ersten 
Systemmode s x zu einer zweiten Systemmode s^ ubergeht, eine 
Folge von gemischten Vorhersagemodellen g L erfafit wird. Die 
erf indungsgemaJJe Vorrichtung kann ferner Einrichtungen zur 
Einstellung von Interpolations- und Segmentierungsparametern, 
Vergleicherschaltungen zur Verarbeitung von Vorhersagef ehlern 
von Vorhersagemodellen, Anzeige- und Signalisierungsein- 
richtungen und Speichereinrichtungen enthalten. Die 
erf indungsgemaiJe Vorrichtung kann ein Monitor flir 
physiologische Daten Oder physikalische oder chemische 
Prozefiparameter sein . 

Mit der Erfindung steht ein Instrument zur Verfiigung, welches 
ein groBes Anwendungspotential in vielen medizinischen, 
wissenschaf tlichen und technischen Bereichen hat. Mit der 
Segmentierung von Signalen und einer entsprechenden 
Identif izierung der zugrundeliegenden Dynamik eroffnen sich 
neue Moglichkeiten der Vorhersage und der Kontrolle auch 
wesentlich nichtstationarer Systeme. 

Anwendungen der Erfindung haben gezeigt, dafl sich 
kontinuierliche Obergange zwischen Systemmoden sicher 
identif izieren lieflen und daft die zugrundeliegenden Dynamiken 
durch die Modelle mit einer Prazision beschreibbar sind, die 
in vielen Fallen eine Vorhersage des Systemverhaltens 

ermoglicht. Mit der Erfindung lassen sich in vielen Fallen • 
nichtstationarer Prozesse Modelle identif izieren, die sich 
auch fiir eine Kontrolle der Prozesse eignen, die ohne eine 
Beriicksichtigung der Nichtstationaritat nicht moglich waren. 
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Ausfuhrungsf ormen und weitere Vorteile der Erfindung werden im 
folgenden unter Bezug auf die beigefiigten Zeichnungen 
beschrieben. Es zeigen: 

Fig. 1: Kurvendars tellungen zur Illustration eines ersten 
Segmentierungsschrittes des er f indungsgemalien Verfahrens; 

Fig. 2: Kurvendarstellungen zur Illustration eines weiteren 
Segmentierungsschrittes des er f indungsgemaften Verfahrens 

Fig. 3: Kurvendarstellungen zur Segment ierung von 
Blutregulierungsdaten nach dem erf indungsgemaiien Verfahren; 
und 

Fig. 4: Kurvendarstellungen zur Segmentierung von EEG-Daten 
nach dem erf indungsgemaJien Verfahren. 

Im folgenden werden zunachst Einzelheiten der Erfindung unter 
Bezug auf die Figuren 1 und 2 und anschliefiend praktische 
Anwendungsbeispiele erlautert. Es ist dem Fachmann 
ersichtlich, dafi die Erfindung nicht auf die 

Anwendungsbeispiele beschrankt ist, sondern entsprechend auch 
in anderen Gebieten angewendet werden kann, wie sie 
beispielhaft weiter unten aufgefuhrt sind. 

(1) Erfassung von Drif tubergangen in nichtstationaren 
Zeitreihen 

Erf indungsgemaJi werden nichtstationare Zeitreihen mit einer 
zweischrittigen Prozedur erfafit, bei der zunachst eine 
geeignete Modellierung und anschliefiend eine sogenannte 
Drif tsegmentierung erfolgt. Die Modellierung ist dazu 
eingerichtet/ in jedem Zeitabschnitt einer vorbestimmten 
Mindestlange fiir jeden Systemparameter ein vorbestimmtes 
Vorhersagemodell fiir eine entsprechende Systemmode zu 
erfassen. Hierzu erfolgt vorzugsweise eine herkdmmliche 
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Schal tsegmentierung, wie sie beispielsweise aus der 
Publikation von K. Pawelzik et al . in "Neural Computation" 
(Bd. 8, 1996, S. 340 ff.) bekannt ist. Die Modellierung kann 
aber auch durch eine andere, in Bezug auf die gewonnene 
Systeminf ormation zur Schaltsegmentierung aquivalente Prozedur 
erfolgen, die an eine konkrete Anwendung z. B. bei bekannten 
reinen Moden oder Randbedingungen angepaiit ist. 

Die Schritte der Schalt- und Drif tsegmentierung werden im 
folgenden naher erlautert. Hinsichtlich der 

Schaltsegmentierung wird hiermit der Inhalt der Publikation 
von K. Pawelzik et al. ausdrticklich vollstandig in die 
vor liegende Beschreibung einbezogen . 

(i) Schritt 1 (Schaltsegmentierung) 

Die Schaltsegmentierung dient der Ermittlung 
charakteristischer Prediktoren, die zur Beschreibung der 
Systemmoden geeignet sind. Die Schaltsegmentierung kann 
entweder an einer Tr ainings-Zeitreihe oder an der zu 
untersuchenden Zeitreihe durchgefiihrt werden. In beiden Fallen 
konnen die ermittelten Vorhersagemodelle oder Prediktoren fur 
weitere, unbekannte Zeitreihen verwendet werden. 

Es wird ein dynamisches System mit einer endlichen Zahl N 
unterschiedlicher Moden betrachtet. Fur die j-te Mode ist ein 

Wert (allg.: Vektor oder Satz) ot^t) eines beobachtbaren 
Systemparameters charakteristisch, der jeweils mit einer 
Funktion f i(t) (i = 1,...,N) aus einem Satz von N Funktionen f 
modelliert werden soli. Es wird die Zeitreihe {x t } = x^(t) der 
Systemvariablen betrachtet und zeitabhangig jeweils nach der 
Funktion f i(t) gesucht, fur die {y t } = y^t) = f i(t) (x j (t) ) eine 
neue Zeitreihe vorher zusagender Punkte y^Ct) darstellt, die in 
Bezug auf die Systemmoden qualitativ die selben Eigenschaf ten 
hat wie {x t } . Durch den Wechsel der Modellf unktion f im 
Zeitverlauf wird die Schaltsegmentierung gefunden, die die 
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Zeitreihe {x t } entsprechend der wechselnden Systemmoden 
unterteilt . 

Die Funktionen f werden als Prediktoren (oder: 
Vorhersagemodelle, Expertenf unktionen) aus einem Satz von 
Netzwerken mit verander lichen Parametern durch ein geeignetes 
Trainingsprogramm abgeleitet, bei dem simultan sowohl die 
Parameter der Netzwerke als auch die Segmentierung ermittelt 
werden. Der Begriff "Netzwerk" wird hier fur alle moglichen 
geeigneten Modellf unktionen verwendet, also vorzugsweise fur 
neuronale Netzwerke, aber auch z. B. ftir Polynome oder lineare 
Funktionsapproximationen . Die optimale Wahl eines neuronalen 
Netzwerkes erfolgt in Abhangigkeit von der spezifischen 
Anwendung . Es werden vorzugsweise Netzwerke mit einer 
schnellen Lernf ahigkeit , wie z. B. sogenannte RBF-Netzwerke 
(Radial Basis Function Network) vom Moody-Darken-Typ 
verwendet . 

Das Training erfolgt unter der Voraussetzung, daii die 
Systemmoden nicht mit jedem Zeitschritt wechseln, sondern eine 
geringere Schaltrate aufweisen, so daii eine Systemmode liber 
mehrere Zeitschritte erhalten bleibt. Die angenommene Grenze 
der Schaltrate bzw. Anzahl der Zeitschritte, iiber die eine 
Systemmode erhalten bleibt, ist zunachst ein freier 
Eingangsparameter und kann anwendungsabhangig in geeigneter 
Weise beispielsweise in Abhangigkeit vorgegebener 
Erf ahrungswerte oder nach einer Parameteranpassungsstrategie 
ausgewahlt werden, Bei der Parameteranpassungsstrategie kann 
vorgesehen sein, fur die Schaltrate einen Ausgangswert 
vorzugeben und mit diesem einen Vorhersagef ehler (siehe unten) 
zu ermitteln. Falls die Schaltrate zu hoch oder zu niedrig 
gewahlt ist, ftihrt eine entsprechende Ober- oder 
Unterspezialisierung zu einem zu hohen Vorhersagef ehler . Im 
weiteren Verlauf der Anpassung kann dann die Schaltrate 
optimiert werden, bis der mittlere Vorhersagef ehler unterhalb 
vorbestimmter Grenzen liegt. 
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Das Training erfolgt durch Maximierung der Wahrscheinlichkeit 
W, da!3 der Satz von Netzwerken die Zeitreihe {x t } erzeugen 
wiirde . Es handelt sich urn ein Training mit konkurrierendem 
Lernen (engl.: "competitive learning"), wie es im einzelnen in 
der Abhandlung "Introduction to the theory of neural 
computation" von J. Hertz et al . (Addison-Wesley Publishing 
Company 1991, insb. Kap. 9: "Unsupervised competitive 
learning") beschrieben wird. Aus dieser Abhandlung ist die 
anwendungsabhangige Implementierung eines solchen Trainings 
ableitbar. Die Trainingsregel des konkurrierenden Lernens auf 
der Grundlage des beim Lernen auftretenden Fehlers ist gemali 
(1) darstellbar. 



Diese Trainingsregel stellt sicher, dafi die Lerngeschwin- 
digkeit (Verbesserung der Parameter) fur die Funktionen f mit 
geringstem Abstand zum Zielwert y am hochsten ist. 

Figur 1 zeigt das Ergebnis der Schaltsegmentierung am Beispiel 
der Analyse einer chaotischen Zeitreihe {x t } mit x t+1 = f (x t ) 
zwischen den vier Moden: 

f x (x) = 4x(l-x) fur x e [0, 1] 
f 2 (x) = f, (f x (x) ) 

f 3 (x) = 2x ftir x g [0, 0.5) bzw. 

f 3 (x) = 2(l-x) fur x e [0,5, 1] 
f 4 (x) = f 3 (f 3 (x) ) 

Ftir die ersten 50 Zeitschritte wird mit einem Startwert 
x 0 = 0.5289 zunachst f, angewendet . Anschliefiend erfolgt ein 
Obergang (Einzelheiten siehe (ii) ) zur Mode f 2 , die nach 
Schritt 100 bis Schritt 150 stationar wird. Entsprechend wird 




e -P(y-fi) 2 



(y-f±) 




-P(y-fi) 2 



(1) 
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ab Schritt 200 bzw. Schritt 300 jeweils fur 50 Schritte die 
Mode f 3 bzw. f 4 eingenommen . Danach erfolgt der Riickiibergang 
auf f x . Figur la zelgt einen Ausschnitt (Schritte 300 bis 450) 
des Zeitverlaufs der Zeitreihe {x t } mit x t+1 = f (x t ) . 

Die Segmentierung der ersten 450 Zeitschritte mit 6 
Prediktoren ^ i = 1,...,6 (RBF-Net zwerke vom 
Moody-Darken-Typ) ist in Figur lb dargestellt. Das Training 
ergibt eine Spezialisierung von vier der Prediktoren (6, 2, 4, 
3) jeweils auf eine der o. a. vier Moden. Die stationaren, 
Bereiche liegen bei den Intervallen [0, 50] und [400, 450] 
(fj, [100, 150] (f 2 ), [200, 250] (f 3 ) und [300, 350] (f 4 ). Die 
ubrigen zwei Prediktoren (3, 5) haben sich auf die 
Ubergangsbereiche zwischen den Moden spezialisiert . Dies zeigt 
den Nachteil der herkommlichen Schal tsegmentierung, bei der im 
Fall von Ubergangen der entsprechende Zeitbereich ohne 
adaquate Beschreibung mehrfach unterteilt wird. 

Anstelle des hier beschriebenen sog. "Hard competition"- 
Trainings, bei dem bei einem Trainingsschritt jeweils nur ein 
Vorhersagemodell optimiert wird ("winner takes all") kann auch 
vorgesehen sein, den Grad der Konkurrenz im Rahmen eines sog. 
"Soft competition"-Trainings zu verandern, wie es im einzelnen 
in der Publikation von K. Pawelzik et al. beschrieben ist. 

(ii) Schritt 2 (Drif tsegmentierung) 

Beim zweiten Schritt werden die Ubergange (sog. Drif ten, 
nicht-abruptes Cbergehen, gleitendes Wechseln) zwischen den 
Systemmoden beriicksichtigt . Mit der Erfindung wurde als 
wichtige Voraussetzung fur die Drif tsegmentierung 
herausgef unden, dafi der Ubergang von einer ersten Systemmode 
direkt in eine zweite Systemmode und nicht iiber eine dritte 
Systemmode erfolgt. Das Drif ten zwischen Systemmoden wird 
somit wie folgt durch eine Uberlagerung von (oder paarweise 
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lineare Interpolation zwischen) genau zwei Moden modelliert. 
Dabei treten gemischte, ggf. abgestufte Zwischenmoden auf, die 
jedoch nicht eigene (reine) Systemmoden sind. 

Es wird ein Satz von P reinen Systemmoden, die jeweils durch 
ein Netzwerk k(s), s e P, reprasentiert werden, und ein Satz 
von M gemischten Systemmoden betrachtet, die jeweils durch 
eine lineare Uberlagerung von zwei Netzwerken i(s) und j (s) , 

s g M, reprasentiert werden. Das Modellnetzwerk g 3 fiir eine 

gegebene Mode s e S, S = P^M wird gemafi Gleichung (2) 
dargestellt . 



Prediktoren, die gemaft der o. a. Schaltsegmentierung ermittelt 
wurden. Dabei ist m eine Einbettungsdimension und x der 
Verzogerungsparameter der Einbettung. Die Einbettungsdimension 
ist die Dimension des Phasenraumes, in dem das System 
betrachtet wird und in dem die Modelle operieren. 

Fiir jede gemischte Systemmode sind zwei Parameter a, b 
zusammen mit zwei Netzwerkindizes i, j charakteristisch ♦ Zur 
Vereinf achung des Rechenauf wandes wird die Zahl der gemischten 
Moden begrenzt. Mit 0 < a(s) < 1 und b(s) = 1 - a(s) wird 
eine endliche Zahl von Werten a(s) definiert. Zur weiteren 
Vereinf achung werden gleiche Abstande zwischen den Werten a(s) 
gemaii Gleichung (3) gewahlt. 

a r = r / (R +1) mit r = 1, . . . , R (3) 

R entspricht der Anzahl der zugelassenen Zwischenmoden und 
wird auch als Auflbsung oder Abstufung der Interpolation 




fiir s g P 
+ b (s) f j(s) (3T t ) fur s e M 



(2) 



In (2) ist x der Vektor (x t , x t _ t , . . . , x t-<m-i>T> der 
Zeitverzogerungskoordinaten der Zeitreihe {x t } und f^ sind 
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zwischen den reinen Moden bezeichnet. Die Auflbsung R kann 
einen beliebigen Wert annehmen, wird jedoch anwendungsabhangig 
zur Erzielung einer optimalen Systembeschreibung ( insbesondere 
bei stark verrauschten Vorgangen) und praktikabler 
Rechenzeiten insbesondere unter Berucksichtigung der oben 
genannten Schaltrate ausreichend niedrig gewahlt. Bei 
praktischen Anwendungen (siehe unten) kann vorgesehen sein, 
daii die Auflbsung R manuell von einem Bediener oder 
automatisch von einem Stellkreis in Abhangigkeit von einem 
vorliegenden Analyseergebnis und einem Vergleich mit einem 
vorgegebenen Schwellwert gewahlt werden. 

Die Gesamtzahl der gemischten Moden betragt bei einer 
gegebenen Auflbsung R zwischen jeweils zwei Netzen 
I m I = R ' N * (N-D/2. Bei dem o. a. Beispiel betragt somit bei 
N = 8 reinen Moden und einer Auflbsung R = 32 und Gesamtzahl 
der gemischten Moden |M| = 896. Zur Bestimmung der Gesamtzahl 
der Systemmoden kommen noch die 8 reinen Moden hinzu. 

Die Drif tsegmentierung umfafit nun die Suche nach einer 
Segmentierung mit den reinen und gemischten Systemmoden (a, b, 
R) , die in Bezug auf den Vorhersagef ehler der Moden der 
gesamten Zeitreihe optimiert ist. Die Prediktoren werden so 
ausgewahlt, daii jedem Element der Zeitreihe einer der Moden 
aus der Gesamtzahl der Systemmoden zugeordnet werden kann. Der 
Vorhersagef ehler ist die Abweichung einer 
Prediktorenvorhersage vom tatsachlichen Element der zu 
untersuchenden Zeitreihe. Fur die zu untersuchende Zeitreihe, 
die nicht mehr unbedingt die Trainings-Zeitreihe ist, mit der 
bei der Schaltsegmentierung die angepaflten Netzwerke oder 
Prediktoren ermittelt wurden, wird fur jeden Zeitschritt mit 
jedem der Prediktoren eine Vorhersage ermittelt, woraus sich 
eine zeitabhangige Matrix der Prediktorenvorhersagen ergibt, 
aus der ein mittlerer Vorhersagef ehler fur willkurlich 
gewahlte Segmentierungen ableitbar ist. Die Segmentierungen 
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mit dem geringsten Vorhersagef ehler ist die gesuchte 
Drif tsegmentierung . 

Die Suche nach der Segmentierung mit dem geringsten 
Vorhersagef ehler kann mit jeder geeigneten Such- oder 
Iterationstechnik erfolgen. Vorzugsweise wird eine dynamische 
Programmiertechnik gewahlt, die Equivalent zum 

Viterbi-Algorithmus fiir HM-Modelle (sog. Hidden Markov Models) 
ist. Einzelheiten hierzu sind beispielsweise in der 
Publikation "A Tutorial on Hidden Markov Models and Selected 
Applications in Speech Recognition" von L. R. Rabiner in 
"Readings in Speech Recognition" (Hrsg. A. Waibel et al . , San 
Mateo, Morgan Kaufmann, 1990, S. 267-296) beschrieben. Die 
Drif tsegmentierung ist im Rahmen der HM-Modelle die 
wahrscheinlichste Modensequenz, die die zu untersuchende 
Zeitreihe erzeiigt haben konnte. Als Nebenbedingung wird dabei 
die Mbglichkeit der Moden-Anderungen iiber die T-Funktion 
(siehe unten) eingeschrankt . 

Das Ziel der Anpassung ist die Angabe einer optimalen Sequenz 
von Netzwerken bzw. linearen Mischungen aus diesen. Eine 
Sequenz ist optimal, wenn die sogenannte Energie- oder 
Cost-Funktion C* der Vorhersage minimiert ist. Die 
Cost-Funktion C* setzt sich aus der Summe der quadratischen 
Fehler der Vorhersage und den Cost-Funktionen der 
Modeniibergange der Sequenz zusammen. Die Ableitung der 
Cost-Funktion C* jeweils zwischen zwei Zeitpunkten t 0 und t max 
erfolgt induktiv, indem zunachst von einer Start-Cost-Funktion 
gemaB Gleichung (4) ausgegangen wird. 



C s (t 0 ) = e s (t 0 ) 



(4) 



wobei 



s s (t) = <x t - g.ilt^)) 2 



(5) 



der quadratischen Fehler der Vorhersage der reinen oder 
gemischten Moden g ist. 
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Fur cien Induktionsschr itt von t - 1 auf t wird die 
Cost-Funktion gemafl Gleichung (6) fur alle s e S berechnet: 

C s (t) = e s (t) + min {C s (t-1) + T(s, s)}, t=t 0 +l, . . . t maK (10) 

Dabei ist T(s, s) die Cost-Funktion des Ubergangs von einer 
Mode S zu einer Mode s . 

Die optimale (minimale) Cost-Funktion C* ist dann: 

C* = min (C s (t max )} (11) .„ 

Die Funktion T entspricht bei den HM-Modellen den 
Obergangswahrscheinlichkeiten und kann anwendungsabhangig 
geeignet gewahlt werden. Es kann beispielsweise vorgesehen 
sein, abrupte Schaltlibergange und gleitendes Driften zwischen 

zwei Netzen zuzulassen und alle anderen Ubergange durch T = oo 
auszuschlieiien . 

Die Drif tsegmentierung ergibt sich aus der ermittelten 
optimalen Sequenz von Netzwerken bzw. linearen Mischungen aus 
diesen, indem die Moden, die C* ergeben, zuriickverf olgt und in 
Abhangigkeit von der Zeit erfailt werden. 

Im Anschlufi an die Drif tsegmentierung kann sich ein 
zusatzlicher Schritt der Reduzierung der Zahl der zur 
Modellierung verwendeten Netzwerke anschlielien, der im 
einzelnen unten beschrieben wird. 

Abschliefiend findet eine Identif izierung der segmentierten 
Moden statt, indem jedem Prediktor oder Vorhersagemodell die 
zugehorige Systemmode zugeordnet wird. Diese Identif izierung 
erf olgt anwendungsabhangig . 
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Das Ergebnis der Drif tsegmentierung im Fall der oben unter 
Bezug auf Figur 1 erlauterten chaotischen Zeitreihe {x t } mit 
vier Moden wird im folgenden unter Bezug auf Figur 2 
beschrieben. Die Drif tsegmentierung umfafit die Suche nach 
einem Verlauf a(t), der einen speziellen Pfad zwischen den 
reinen Moden liefert, fur den der Vorhersagef ehler der 
gesamten Zeitreihe optimiert ist. 

Jeweils nach den ersten 50 Zeitschritte mit der Mode gemafi f 1 
erfolgt fur 50 Schritte ein zeitlinearer Ubergang in die Mode 
gemafi f 2 . Der Obergang ist ein zeitabhangiges Drif ten gemafi den 
Gleichungen (12) : 



Entsprechende Obergange finden nach dem 150., 250. und 350. 
Schritt fur jeweils 50 Schritte statt. 

In Figur 2 ist die Besetzung der jeweiligen Moden entsprechend 
den ermittelten Netzwerken in Abhangigkeit von der Zeit 
(Zeitschritte [1200, 2400]) dargestellt. Aus 
Obersichtlichkeitsgrunden sind die Obergangs- oder 
Drif tbereiche entsprechend ihren Zeitgrenzen und Ausgangs- 
bzw. Endmoden in Rahmen dargestellt, in denen der jeweilige 
Driftverlauf zwischen den Moden gepunktet ist. Figur 2a zeigt 
bei einer Auflosung R = 32 (siehe Gleichung (3)) Obergange wie 
z. B. bei den Zeitschritten 1350 bis 1400 zwischen den Netzen 
2 und 4. Die Obergange sind wie gemafi Gleichung (12) erwartet 
linear. Bei einer geringeren Auflosung von R = 3 ergibt sich 
die in Figur 2b dargestellte Segmentierung . Abweichend vom 
linearen Driftverlauf sind die gepunkteten Obergange 
stuf enformig. Dennoch ist die Darstellung auch bei der 
geringeren Auflosung eine adequate Beschreibung des 



f("x t ) = (1 - attllf^^J + a(t)f 2 ("x t ) 



mit 



(12) 




= 100 
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dynamischen Verhaltens des Systems, wie der Vergleich der 
zeitlichen Lage der Moden und des Driftens zeigt. 

(2) Anwendungsbeispiele ftir die Erfassung von Drif ttibergangen 

(i) Blutzellenregulation im menschlichen Korper 

Die Blutzellenregulation im menschlichen Korper stellt ein 
hochdimensionales chaotisches System dar, das durch die 
Mackey-Glass-Verzogerungs-Dif f erentialgleichung (13) (siehe 
auch in der o. a. Abhandlung von J. Hertz et al.) beschrieben 
werden kann. 

dx(t)/dt = -O.lx(t) + (0.2x (t-td) / (1+x (t-t d ) 10 ): (13) 

Erfindungsgemali konnen Zeitreihen von physiologischen 
Parametern, die ftir die Menge der roten Blutkorper 
charakteristisch sind, anwendungsabhangig segmentiert werden. 
Die Funktionsf ahigkeit der Segment ierung wird im folgenden 
beispielhaf t erlautert . 

Bei Vorgabe von zwei Moden A und B, die sich jeweils durch 
verschiedene Verzogerungsparameter t d = 17 bzw. t d = 23 
unterscheiden, erfolgt mit einer Abtast-Zeitschrittgrofte x = 6 
nach 100 Schritten zunachst ein Obergang von A nach B. Der 
Obergang dauert 100 Schritte und ist eine Uberlagerung der 
Gleichung (13) mit den beiden Verzogerungsparameter t d wahrend 
der Integration von Gleichung (13) . Die Uberlagerung wird mit 
einem exponentiellen Drif tparameter a (siehe Gleichung (2) ) 
gemaft Gleichung (14) erzeugt. 

a(t) = exp(-4t/100) , t - 1,..., 100 (14) 

In der Folge wiederholen sich alle 100 Schritte stationare 
Moden A oder B bzw. die jeweiligen Obergange. Fur jeden 
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Ruck-Obergang nach einem Drif tlibergang wird ein schaltartiges 
Umspringen angenommen . Figur 3a zeigt die entsprechende 
Zeitreihe liber 300 Schritte. Die Dri f tsegmentierung mit sechs 
Prediktoren auf der Grundlage von RBF-Net zwerken mit jeweils 
4 0 Basisf unktionen, einem Einbettungsparameter m = 6 und dem 

Verzogerungsparameter x = 1 (siehe Gleichung (2)) ergibt das 
in Figur 3b gezeigte Bild. Es zeigt sich die erwartete 
Segmentierung der Zeitreihe in stationare Moden und 
Drif tiibergange . 

Allerdings zeigt sich, daft sich jeweils 2 Netzwerke auf eine 
Mode spezialisiert haben (2, 3 => Mode A, 5, 6 => Mode B) . In 
einer solchen Situation kann es er f indungsgemafl vorgesehen 
sein, den zusatzlichen Schritt der Reduzierung der Zahl der 
zur Modeilierung verwendeten Netzwerke vorzusehen. 

Der Reduzierungschritt umfafit eine sequentielle Reduzierung 
der Netzwerkzahl, jeweils verbunden mit einer Bestimmung des 
mittleren Vorhersagef ehlers . Die Reduzierung (Abzug von 
Redundanznetzwerken) wird beendet, falls eine weitere 
Verringerung der Netzwerkzahl eine signifikante Erhdhung des 
Vorhersagef ehlers bedeutet . Figur 3c zeigt das Ergebnis einer 
solchen Reduzierung. Der mittlere quadratische Fehler RMSE 
bleibt bei der Reduzierung urn ein, zwei, drei und vier 
Netzwerke jeweils konstant, wohingegen ein starker Anstieg fur 
den Fall auftritt, dafi. nur noch mit einem Netzwerk modelliert 
wird. Dies bedeutet, dafl das System optimal mit einer Zahl von 
Netzwerken modelliert wird, die gleich der Gesamtzahl von 
betrachteten Netzwerken, vermindert um die Zahl von 
Redundanznetzwerken ist . 

Die adequaten Modellnetze werden dann erhalten, indem der 
RMSE-Wert fur jede Netz-Kombination mit reduzierter Netzzahl 
berechnet wird. Die Netz-Kombination mit dem geringsten 
RMSE-Wert umfafit die gesuchten Modellnetzwerke oder 
Prediktoren. Figur 3d zeigt die Drif tsegmentierung nach dem 
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Reduzierungsschritt . Demnach beschreiben die iibrig gebliebenen 
Prediktoren 2 und 5 das System vollstandig. 

( ii ) Schlaf datener f as sung 

Eine weitere Anwendung der Erfindung liegt im Bereich der 
Analyse physiologischer Daten, die fur den Ablauf von Schlaf- 
und Wachmoden von Lebewesen charakter istisch sind. Als 
Grundlage fur sich anschliefiende Prozeduren zur Erfassung von 
Schlaf storungen konnen Zeitreihen z. B. von EEG-Daten 
segmentiert werden . 

Figur 4a zeigt im Vergleich die Ergebnisse einer her kommlichen 
Schaltsegmentierung (oben) , einer Drif tsegmentierung (Mitte) 
und einer "manuellen" Segmentierung (unten) eines 
medizinischen Fachmannes (Schlaf forscher ) auf der Grundlage 
von Erf ahrungswerten am Beispiel eines Nachmittagsschlaf es 
eines gesunden Menschen. Die Schalt- und Drif tsegmentierungen 
erfolgen mit acht Netzwerken (netl . .. net8) an 
. Einkanal-EEG-Daten x(t) (Figur 4b). In Figur 4a sind wie bei 
Figur 2 aus Obersichtlichkeitsgrunden Rahmen gezogen, die bei 
den Driftmoden verdeutlichen, zwischen welchen Netzen 
interpoliert wird. Die gepunktete Linie im Inneren der Rahmen 
zeigt den jeweils tatsachlichen Verlauf . Die manuelle 
Segmentierung basiert auf der Beobachtung von physiologischen 
Signalen (z. B. EEG, EOG, ECG, Puis, Blutdruck, Atmung, 
Augenbewegungen) . Die Moden Wl, W2 bezeichnen zwei Wachmoden 
mit offenen bzw. geschlossenen Augen und die Moden SI, S2 
jeweils Schlaf zustande . "n.a." und "art." beziehen sich auf 
nicht in Betracht gezogene Zustande bzw. Artefakte. 

Die Schaltsegmentierung zeigt ein vergleichsweise 
undif f erenziertes Bild, das nur grob mit den ubrigen 
Beobachtungen konsistent ist. So tritt beispielsweise ein 

Einschlafvorgang in alien drei Fallen bei t « 7000 auf. Die 
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Drif tsegmentierung ergibt jedoch mehrere Drif tubergange, die 
zusatzliche Einzelheiten des Schlaf verhaltens darstellen. Der 
"manuell" beobachtete Schlafbeginn zur Zeit t « 4000 wird 
durch einen exponentiellen Drif tiibergang von Netz net7 
(Wachmoden-Prediktor ) zum Netz net4 ( Schlafmoden-Prediktor ) 
reprasentiert . Das Aufwachen setzt bei t ~ 9000 durch ein 
leichtes Ruckdriften zum Netz net7 ein, das bis zum Erreichen 

des "manuell" ermittelten Auf wachpunktes t « 9500 gehalten 
wird. In dieser Situation erfolgt eine plotzliche Anderung des 
Wichtungsf aktors, so daii das Netz net7 eine groBere Wichtung 

erhalt. Nach t » 9800 (Augen of fen) besteht eine Mischung der 
beiden Wachmoden-Prediktoren net7 und net2 . 

(iii) Weitere Anwendungen und Vorteile 

Figur 4a zeigt, daft mit dem erf indungsgemafien Verfahren 
automat isch detaillierte Segment ierungen erzielt werden 
konnen, die bi slang nur durch die Beobachtung komplexer 
Merkmalsbilder auf der Grundlage breiter Erfahrungen und 
Intuitionen zuganglich waren. Dieser Vorteil ist nicht nur in 
der Medizin, sondern auch in anderen Gebieten nutzbar, in 
denen grofie Datenmengen bei der Beschreibung komplexer 
dynamischer Systeme anf alien. Derartige Gebiete sind die 
physikalische, chemische und/oder biologische 
Verf ahrenstechnik, die Geologie/ Meteorologie, Klimatologie, 
die Spracherf assung u. dgl . . 

Erf indungsgemaBe Verfahren haben die folgenden Vorteile. Das 
betrachtete System kann hochdimensional sein (10 oder mehr 
Dimensionen) . Die Erfindung erlaubt eine Reduzierung der 
Komplexitat eines solchen Systems durch Betrachtung von 
niedrigerdimensionalen Moden und von wechselnden Ubergangen 
zwischen diesen. Der Einsatz von Vorhersagemodellen fur die 



20 



WO 99/14686 



PCT/EP98/05793 



Segmentierung ist invariant gegen Anderungen der Amplitude 
erfaflter Signale. 

Die Anwendung der Erfindung fiir die Vorhersage oder Kontrolle 
eines Systems erfolgt derart, daft zunachst wie oben 
beschrieben aus der vergangenen Beobachtung und der Kenntnis 
der aktuellen Moden der IST-Zustand des Systems erfaftt wird, 
der ggf . eine Mischung entsprechend dem Ergebnis der 
Drif tsegmentierung darstellt. Der IST-Zustand entspricht einem 
dynamischen System f. Die Vorhersage bedeutet, dafi das System 
f auf den momentanen Zustand x angewendet wird und sich daraus 
die Vorhersage fiir den unmittelbar folgenden Zustand y ergibt. 
Die Kontrolle bedeutet, dafi aus dem IST-Zustand die Abweichung 
von einem SOLL-Zustand ermittelt und aus der Abweichung eine 
geeignete Regelstrategie abgeleitet wird. 

Der Vorteil der Vorhersage und Kontrolle besteht darin, dali in 
komplexen Systemen (z. B. bei der Erfassung chemischer 
Reaktionen in einem Reaktor) , die ggf. die Messung lediglich 
weniger Mefigroiien erlauben, die fiir sich allein aufgrund von 
Mehrdeutigkei ten oder wegen systemimmanenten Verzogerungen 
keine direkten Ruckschliisse auf den Systemzustand und ggf. 
bestehende gemischte Zustande erlauben, dennoch detaillierte 
Inf ormationen liber das System abgeleitet werden konnen. So 
kann bei dem Beispiel mit einer chemischen Reaktion aus der 
erf indungsgemaJien Erfassung bespielsweise der makroskopischer 
thermodynamischer Zustandsgrofien eine optimale Regelstrategie 
abgeleitet werden, die die Dosierung bestimmter 
Reaktionspartner umf alit . 
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Patentanspriiche 



1. Verfahren zur Erfassung der Moden eines dynamischen 
Systems mit einer Vielzahl von Moden s if die jeweils einen Satz 
a(t) charakteristischer Systemparameter besitzen, wobei eine 
Zeitreihe mindestens einer Systemvariablen x(t) einer 
Modellierung unterzogen wird, die dazu eingerichtet ist, in 
jedem Zeitabschnitt einer vorbestimmten Mindestlange fur jede 
Systemvariable x(t) ein vorbestimmtes Vorhersagemodell f, fur 
eine entsprechende Systemmode s i zu erfassen, 

dadurch gekennzeichnet , dafi 

nach der Modellierung der Zeitreihe eine Drif tsegmentierung 
erfolgt, bei der in jedem Zeitabschnitt, in dem das System von 
einer ersten Systemmode s t zu einer zweiten Systemmode s^ 
iibergeht, eine Folge von gemischten Vorhersagemodellen g L 
erfafit wird, die durch eine lineare, paarweise Uberlagerung 
der Vorhersagemodelle f i#j der zwei Systemmoden s i#j gegeben ist. 

2. Verfahren gemaii Anspruch 1, bei dem die Modellierung eine 
Schaltsegmentierung ist . 

3. Verfahren gemaB Anspruch 2, bei dem die 
Schaltsegmentierung durch Simulierung einer 

Trainings-Zeitreihe des Systems oder der zu untersuchenden 
Zeitreihe mit mehreren, untereinander in Konkurrenz stehenden 
Vorhersagemodellen erf olgt . 

4. Verfahren gemafJ Anspruch 3, bei dem die Vorhersagemodelle 
durch neuronale Netzwerke oder andere Modelle zur Schatzung 
von Funktionen gebildet werden, die jeweils fur eine Mode s 
charakteristisch sind und nach vorbestimmten Trainingsregeln 
urn die Beschreibung der einzelnen Elemente der Zeitreihe 
konkurrieren . 
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5 Verfahren gemafl einem der Anspriiche 1 bis 4, bei dem die 
Folge gemischter Systemmoden g i aus den Vorhersagemodelle f i#j 
und Interpolationsparametern a, b gemafi 
gi = a(s)fi (S )(x) + b (s) f j(s) (x) ermittelt wird. 

6. Verfahren gemaii Anspruch 5, bei dem die 

Interpolationsparameter gemafl 0 < a(s) < 1 und b(s) = 1 - a(s) 
gewahlt sind. 

7. Verfahren gemafi Anspruch 6, bei dem die Werte a(s) auf 
eine bestimmte Auf losungszahl R beschrankt und/oder 
aquidistant sind. 

8. Verfahren gemafi einem der vorhergehenden Anspriiche, bei 
dem die Folge gemischter Vorhersagemodelle q i erfalit wird, 
indem fur jeden Zeitschritt mit jedem der moglichen 
Vorhersagemodelle jeweils eine Vorhersage ermittelt wird, 
woraus sich eine zeitabhangige Vorhersagematrix ergibt, aus 
der ein mittlerer Vorhersagef ehler fur willkurlich gewahlte 
Segmentierungen ableitbar ist, wobei die gesuchte Folge 
gemischter Vorhersagemodelle q i die Segmentierung mit dem 
geringsten Vorhersagef ehler bzw. der maximalen 
Wahrscheinlichkeit ist . 

9. Verfahren gemafi Anspruch 8, bei dem die Suche nach der 
Segmentierung mit dem geringsten Vorhersagef ehler mit einer 
dynamische Programmiertechnik erfolgt, die Equivalent zum 
Viterbi-Algorithmus fur Hidden-Markov-Modelle ist, wobei eine 
optimale Sequenz von Vorhersagemodellen unter Verwendung einer 
minimierten Cost-Funktion C* der Vorhersage ermittelt wird und 
die Segmentierung induktiv aus der Sequenz von 
Vorhersagemodelle abgeleitet wird. 

10. Verfahren gemafi einem der vorhergehenden Anspriiche, bei 
dem nach der Drif tsegmentierung ein zusatzlicher Schritt der 
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Reduzierung der Zahl der zur Modellierung verwendeten 
Vorhersagemodelle erfolgt, bei dem die Zahl der 
Vorhersagemodelle sequent iell jeweils verbunden mit einer 
Bestimmung des mittleren Vorhersagef ehlers so weit reduziert 
wird, dafi eine weitere Verringerung der Zahl der 
Vorhersagemodelle eine Erhohung des Vorhersagef ehlers 
bedeutet . 

11. Verfahren gemafi einem der vorhergehenden Anspruche, bei 
dem die Zeitreihe mindestens einer der Systemvariablen x(t) 
eine Zeitreihe von physiologischen Parametern umfafit, die mit 
der Mackey-Glass-Verzogerungs-Dif f erentialgleichung 

dx(t)/dt = -O.lx(t) + (0.2x(t-t d ) / (1+x (t-td) 10 ) 
beschrieben wird. 

12. Verfahren gemail einem der Anspruche 1 bis 11, bei dem die 
Zeitreihe mindestens einer der Systemvariablen x(t) eine 
Zeitreihe von physiologischen Parametern umfafit, die fur den 
Ablauf von Schlaf- und Wachmoden charakteristisch sind. 

13. Verfahren gemaii Anspruch 12, bei dem die physiologischen 
Parameter EEG-Signale umfassen. 

14. Verfahren gemafi einem der Anspruche 1 bis 10, bei dem die 
Zeitreihe mindestens einer der Systemvariablen x(t) eine 
Zeitreihe von Sprachsignalen umfafit. 
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